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I. OVERVIEW 


Phase 1 of this Investigation has several objectives: (1) to test a new 
experimental climatological model /variable termed the sponge for potential 
blogeographical , ecological, and climatological applications (the sponge 
Is a measure of moisture availability based on dally temperature maxima 
and minima, and precipitation); (2) to Investigate the feasibility of 
utilizing NOAA/AVHRR meteorological satellite data for vegetation classic 
flcatlon; and (3) to Initiate a vegetation gradient model that utilizes 
climatological (I.e., sponge), biological, and NOAA data that Is ultimately 
applicable to global vegetation stratification and monitoring. 

To accomplish the Initial objective, mean monthly and annual sponge values 
are calculated for 75 Texas locations for the “normal" period of 1941-70. 
Similar values are also compuled for approximately 25 stations along an 
east-west transect across Texas for 1979 and 1980. Results suggest that as 
a generalized climatic index, sponge's simplicity and sensitivity make It 
particularly appropriate for trans- regional blogeographic studies. 

The latter two objectives were approached by acquiring vegetation, clima- 
tological (sponge), and AVHRR pixel data (channels 1 and 2) for 12 locations 
alot.g the east-west Texas gradient. The normalized difference (NO) values 
for the AVHRR data when plotted against the vegetation characteristics (bio- 
mass, net productivity , leaf area) and the sponge values suggest that a 
multivariate gradient model Incorporating AVHRR and sponge data may Indeed 
be useful In global vegetation analysis. 
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II. THE SPONGE; A NEW E r C-^CLIMATOLOGICAL VARIABLE 

The planetary distribution >f types of natural vegetation Is largely a 
function of climate, most especially of spatial variations In energy and 
moisture budgets. As a general rule, climate Is recognized as the pre- 
eminent control of natural vegetation at subcontinental -to global-scales. 

For small areas, geologic, pedologic and other local factors may dominate 
for long, even Indefinite, periods of time. 

Natural scientists such as ecologists and biogeographers have attempted to 
utilize climatic measures and Indices In surveys, stratifications and 
classifications of natural vegetation. Typically, It Is assumed that most 
major ecoreglon and native vegetation-region "boundaries" actually represent 
climatic discontinuities or "breaks". This asscmiptlon, while not entirely 
valid, probably Is reasonably accurate. 

A key problem for such scholars nas been the determination of a simple but 
accurate climatic "Index" (or Indices) which would reflect the primary 
spatial variations In moisture and energy balances and, hence, could be 
applied to the classification of native vegetation. Climate Itself, a 
complexly synergistic synthesis of many variables, does not readily succumb 
to quantitative classification. [For example, what Is the "real" boundary 
between a desert and a (semi-arid) steppe? In fact, of course, there Is no 
abrupt statistical limit but rather a gradual transition from one type Into 
another.] 

Identifying and quantifying climatic factors which explain the distribution 
of natural vegetation Is even more challenging. Temperature and precipitation 
by themselves are poor descriptions of climate and, hence, explainers of 
vegetation distributions (Mather and Yoshloka, 1968). 

It Is therefore essential to select and/or develop climatic Indices which 
Influence vegetation growth and development. Most recent approaches to this 
problem emphasize the Importance of moisture availability (e.g., surpluses 
versus deficits) and, more specifically, evapotranspi ration. Unfortunately, 
evapotranspi ration Is measured at very few places, and even evaporation 
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Itself Is not widely monitored. Numerous models have been developed which 
estimate evaporation from maasurenents of air temperature, average wind 
speed, and net radiation (Penman, 1948; Jensen et a1 , 1970), but their 
usefulness Is constrained by the sparcity of stations which record solar 
radiation and v.1nd speed. 

Simpler evapotranspi ration schemes which require only air temperature and 
precipitation -- both commonly measured around the world -- have been de- 
vised (e.g., Thornthwalte, 1948; Griffiths, 1964; Moe, 1965; Trenchard, 
1976). Thornthwalte's classification made use of mean ninthly temperature 
and precipitation values to generate a moisture Index. Because of Its 
relative simplicity and accuracy, Thornthwalte's approach has been widely 
adopted (e.g.. It Is used to calculate the USDA's Crop Moisture Index). 
However, Its use of average monthly temperatures somewhat limits Its 
sensitivity to variations In continental Ity and altitude. 

An alternate method of relating climate to vegetation Is that of multi- 
variate discriminant analysis of climatic variables to determine their 
relative Influence In a particular ecoreglon (e.g., Blogeocllmatic Units 
of Vancouver Island , Klinka and Nuszdorfer, 1979). While very accurate 
for detailed, site-specific studies, the resulting multiple regression 
equations tend to be (a) cumbersome and lengthy, and (b) less applicable 
to regional and global-scale vegetation classifications. 

Thus climatologists, geographers, and ecologists Interested In large-area 
comparisons have found themselves forced to choose between an approach 
which stratifies climate somewhat too broadly (e.g., Thornthwalte's) and 
another which "hides the forest In the trees", viz., too imich emphasis 
on detail (e.g., Klinka and Nuszdorfer). 

Recently, Trenchard and Artley (1981) developed a new climatological/ 
meteorological variable whose simple form, minimal data requirements and 
accuracy make It an Ideal candidate for application to meso- and macro- 
scale blogeographical , agrocllmatological , and ecological Investigations. 
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This hypothetical ntdlun Is ttrmad the sponge (see Figure 1). Sponge's 
rationale Is suamarlaed as follows (Trenchard and Artley): 

We desired a simple moisture variable with a sound 
physical basis that used coamon meteorological variables, 
was suitable over a broad range of climates, and appli- 
cable to a single station. The result was named sponge. 

Sponge Is described as a simple mediwn with 8 Inches 
of water holding capacity which Is Initialized half-full 
of water on 1 January.* Each day, In accordance with 
the hydrologic cycle, water Is added to the medium from 
precipitation and lost through evaporation. Precipitation 
(both liquid and frozen) Is added at the full amount until 
the layer Is saturated. It Is this sponge like behavior 
which gives the variable Its name. Any additional pre- 
cipitation Is assumed to be lost as run-off or drainage. 
Evaporation occurs at a fraction of the Class A Pan rate, 
the exact proportion being the ratio of the current 
contents to the total capacity of the sponge. Either 
actual or estimated evaporation pan values may be used. 

The dally contents of the sponge are defined as; 

S “ ^-1 ^ '’l ■ ^^1 *s^.t/cap) 


Where: 

» Sponge contents on day 1, in Inches. 

P^ * Precipitation on day 1, In Inches. 

E « Actual or estimated pan evaporation In Inches on day 1. 
CAP > Sponge capacity In Inches 
and 0 ^ S| ^ CAP 

When evaporation pan measurements are not available, they 
may be estimated with a divisor of 30 days to convert the 
evaporation function to a dally value. 

Si - S^_^ + P^ - EP(TX^, TN^)*S^_^/CAP*30 


♦Alternatively, the final value of the previous year may 
be used as an Initial value, and the capacity may be varied 
for a particular region. 
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Uher«: 

EP " Pan tvaporatlon function 
TX^ ■ Maximum t^ptraturc on day 1. 

TN^ - Minimum temperature on day 1. 

Because of Its simple data requirements (dally pre- 
cipitation and evaporation estimated from maximum and 
minimum temperatures), the sponge can be calculated at 
any temperature precipitation observation station. 

Long-term (1941-1970) sponge "normals'* (average values) were recently cal- 
culated for all first-order meteorological stations In the countermlnous 
USA (Trenchard, 1981). Some representative values are presented In Table 1. 
It Is Interesting to note that mean annual sponge values were found to range 
between Yuma, Arizona (0.20") and Nt. Washington, New Hampshire (8.00", or 
absolute sponge capacity). 

Sponge Is concluded to present a meaningful measure of areal "environmental 
moistness". As a generalized climatic Index, sponge's simplicity and 
sensitivity make It particularly appropriate for transreglonal blogeograph- 
Ical studies (e.g., a large-area vegetation classification). 
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TABLE I. SELECTED USA 8-INCH SPONGE NORMALS: 1941-1970 
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Calculations by M. Trenchard (Lockheed), NASA-JSC, Houston, TX., July, 1981. 



III. THE SPONGE VARIABLE AS APPLIED TO MOISTURE AND VEGETATION 
GRADIENTS IN TEXAS 


There exists In Texas the most pronounced continuous, non-orographic, 
Intrastate climatological gradient found anywhere In the United States. 

At least four distinct, first-order climatic types occur within the state 
(humid subtropical, tropical steppe, tropical desert, and mid-latitude 
steppe), with many more important subtypes (e.g., subtropical subhumid). 

In particular, there is an extraordinarily steep east-west moisture 
gradient, ranging from very himid in southeastern Texas (mean annual 
precipitation >50") to true desert in far western Texas (average yearly 
precipitation < 8.0"). This gradient strongly influences ecological 
patterns, and virtually controls the regional distribution of natural 
vegetation. 

Texas is, then, an excellent natural "laboratory" to test the responsive- 
ness and usefulness of the sponge variable (e.g., with respect to the 
classification of natural vegetation utilizing satellite data). With this 
in mind, mean sponge values were calculated for various Texas locations in 
order to address these questions: 

1. Does use of the sponge portray the distribution of climates 
(especially moisture regions) in Texas better than, say, 
precipitation alone? 

2. If sponge accurately reflects the climates of Texas, can these 
values be meaningfully correlated with vegetation-index ("green- 
ness") values as measured from space by NOAA meteorological 
satellites (see the discussion of these indices later in this 
report)? 

3. If the answer to question (2) is affirmative, can a combination 
of sponge and satellite-derived greenness indices be used to 
classify the major natural vegetation regions/types of the 
state? If so, It might well prove feasible to utilize this 
methodology for other large-area and even global -scale vegetation 
surveys and classifications. 
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Climatic Strata And Gradients In Texas Using The Sponge; Long-Term “Normals". 


In order to assess sponge's potential usefulness as a climatic Index In Texas, 
long-term annual and monthly sponje "normal" (1941-1970 mean) values were com- 
puted for 75 locations widely distributed throughout the state (see Map 1). 

Sponge values were obtained by utilizing (1) the formula presented In section 
1 of this report, and (2) mean monthly temperature maxima and minima, and pre- 
cipitation, as compiled by the U.S. National Weather Service (NOAA)*. 

The results of these computations are Illustrated In Maps 2-fi, while Tables 
AI (2) - (f) (pages 48-53) present average normal (1941-70) monthly and annual precip- 
itation and sponge values for each of the 75 locations. They are largely 
self-evident, but several of the more Intriguing aspects should be briefly 
addressed. 

Mean annual sponge values are greatest In the southeast (e.g., 5.2" at 
Beaumont) and decrease continuously to lows In the westernmost quadrant 
(0.49") at Presidio In the Chihuahuan Desert). This Is virtually Identical 
to the pattern of average annual precipitation. However, sponge appears to 
more accurately portray (1) seasonal moisture changes across the state, and 
(2) the magnitude of differences in the relative moistness of the varous 
parts of Texas than does either precipitation or potential evapotranspiratlon 
(estimated by Thornthwalte's method; see Table 3). 

With respect to season. Figures Al - A4 (pages 55-58) and Tables AI (a) - (f) and 
All (page 54) clearly Indicate that, for vegetative activity, winter and spring 
are the wet seasons In East Texas while suiraner and fall are the moist periods In 
West Texas — quite unlike the seasonal distribution of precipitation alone, 
which Is greatest In the summer throughout the state. (Thornthwalte's "Index 
of Moisture" would also reveal this aspect but, because of its reliance on 
mean monthly temperatures, with less spatial sensitivity than sponge; Carter 
and Mather, 1966). 


*Da11y data were simulated from monthly mean temperature maxima, minima , and 
precipitation for each station for 1941-70 using a series of harmonic trans- 
formations. For 1979 and 1980, actual dally data were used. 
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MAP 1. 8- INCH SPONGE NORMALS FOR TEXAS 
(PERIOD 1941-70, EXCEPT AS NOTED IN TABLES AI (a) -(f)) 







HAP 2. NORMAL ANNUAL AVERAGE OF 8-INCH SPONGE - STATE OF TEXAS 
(PERiro 1941-70, EXCEPT AS NOTED IN TABLES AI (a) - (f)) 




















MAP 6. NORMAL OCTOBER AVERA6E 8-lNCH SPONGE - STATE OF TEXAS 
(PERIOD 1941-70, EXCEPT AS NOTED IN TABLES AI (a) - (f)) 






TABLE II. COMPARISON OF NORMAL AVERAGE SPONGE. POTENTIAL EVAPOTRANSPIRATION^ 
ACTUAL EVAPOTRANSPIRATION^ AND PRECIPITATION AT FOUR TEXAS STATION^ 


Station 

Period 

Precipitation 

El Paso/La Tuna 

Annual 

8.06 

El Paso/La Tuna 

January 

0.41 

El Paso/La Tuna 

April 

0.15 

El Paso/La Tuna 

July 

1.61 

El Paso/La Tuna 

October 

0.70 

Me Gamy 

Annual 

12.75 

NeCatney 

January 

0.64 

McCamey 

April 

0.77 

Me Carney 

July 

1.64 

McCamey 

October 

1.39 

Temple 

Annual 

33.87 

Temple 

January 

2.35 

Teny>le 

April 

3.67 

Temple 

July 

1.96 

Temple 

October 

2.73 

Port Arthur/Beaumont 

Annual 

54.77 

Port Arthur/Beaumont 

January 

4.57 

Port Arthur/Beaumont 

April 

4.43 

Port Arthur/BeaufflOnt 

July 

5.71 

Port Arthur/Beaumont 

October 

3.19 


Potential 

Evapotranspi ration 

Actual 

Evapotranspi ration 

Sponge 

38.40 

8.72 

0.59 

0.4f 

0.44 

0.71 

2.72 

0.32 

0.35 

6.96 

1.76 

0.66 

2.60 

0.84 

0.73 

42.72 

14.28 

0.94 

0.32 

0.32 

0.98 

3.24 

0.88 

0.61 

7.48 

1.60 

0.99 

3.04 

1.20 

1.14 

41.20 

33.68 

3.03 

0.44 

0.44 

4.03 

2.92 

2.92 

3.50 

7.40 

5: 20 

2.03 

2.92 

2.92 

2.24 

43.44 

42.44 

5.20 

0.76 

0.76 

7.00 

3.12 

3.12 

5.33 

7.20 

7.08 

4.30 

3.40 

3.36 

3.85 


V. E. and A. E. values extracted from: Average Climatic Water Balance pate of t he Continents. 
Part VII, United States . 1964, C. W. Thornthwaite Assoc Lacoratory of Climatology, in 
dim. Vol. 1) . Wo. 3. Centerton. N. J. 

2 

All values in Inches. 
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Regarding differences In the absolute magnitude of available moisture from 
place to place (I.e., hoM much wetter Is site "x” than site "Y**?), sponge 
also proves highly effective. For example, note that Beaumont's mean pre- 
cipitation Is 6.4 times that of Presidio annually and, In April, Beaumont's 
average rainfall Is 21 times that of Presidio. Sponge shows that the 
moisture gradient between these locations Is actually much steeper : Beau- 

mont's mean annual sponge value Is 10.5 times greater than that of Presidio 
and, In April, Presidio's average sponge value of 0.15" Is only 3* of 
Beaumont (5.33"), a difference of 35. 5X. In other words. West Texas Is 
nearly twice as dry -- compared with the humid eastern part of the state -- 
as precipitation averages alone would suggest. Considering that moisture 
availability Is the primary limiting factor with respect to ecoreglons and 
natural vegetation communities In Texas, It may be concluded that the sponge 
variable is an effective tool for analyzing climate-vegetation relationships. 

Recent Sponge Conditions In Texas; 1979 and 1980 . Mean monthly and annual 
sponge values for the 75 test stations, as well as an additional number of 
locations along the "Texas Transect" (see next section) were calculated for 
1979 and 1980 to assess sponge's responsiveness to Inter-annual moisture 
variability. (Only briefly examined here, these variations, and their 
associations with satellite-measured vegetative index values, will be more 
Intensively studied In a later phase of this research effort.) Refer to 
Tables AI and All and Figures A1 - A4 (pages 48-58). 

It Is evident that 1979 was wetter than normal In East Texas (e.g., 1979 
Ycp at Liberty = 5.31" compared with normal annual of 4.41"), while It 
was dry In central Texas (e.g., at Brady, In 1979 X^p ■ 2.09"; 

Yspnormal = 1.96") and near-normal In West Texas (e.g., at Balmorhea, 

^Sp * - 0.81"). 
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Intrastate moisture conditions In 1980 were quite different than those of 
1979. East Texas was unusually ^ In 1980: For example, Huntsville's mean 
annual sponge was 2.67", only 6U of the long-term normal value. The differ- 
ence was especially pronounced In mid-simmer, when this area experienced 
drought conditlor.s (see Figure A3, page 57). 

By contrast, 1980 was a relatively moist year In West Texas: Pecos and 
Balmorhea, for Instance, had annual sponge values nearly 100X above their 
30-year normals (see Table All, page 54). 

Preliminary Assessment of Sponge 

Based on these early results. It may be concluded that the sponge Is a 
useful new climatic variable for purposes of Identifying and Interpreting 
trans-reglonal moisture (and, therefore, ecological) gradients and strata. 

In fact. It may well prove to be, on bala ice, the best such measure yet 
devised for practical 1arge-area analysis. Accordingly, sponge Is utilized 
In the following sections of this report as a generalized climatic Index, 
one which Is correlated with vegetative Indices derived from NOAA meteoro- 
logical satellite Imagery, as part of a gradient study of natural vegetation 
along a hypothetical east-west "Texas Transect". 
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IV. THE VEGETATION 6RADIEHT IfflLIZINe HOM SATELLITE IMAGERY 


Historically, two broadly concalvad rattarch mtthodi havt avolvad to allow 
stratification and abstraction of pi ant cMmunI t1 as . classification and 
gradient analysis (Kessell, 1979). Classification Involves grouping samples 
together on the basis of shared characteristics Into an abstract dess of 
plant communities. Such a grouping of ccMwiunItles by any definition of 
shared characteristics Is referred to as a community- type (Mhittaker, 1975). 

The secono mrthod, gradient analysis, deals not with discontinuous classes 
but with continuity and gradient relationships. When the arrangement Is 
along a predetermined environmental gradient. I.e. moisture, the method Is 
termed direct gradient analysis . Indirect gradient analysis Is the arrange- 
ment of samples along abstract axes that may or may not correspond to envir- 
onmental gradients. The process of arranging samples along one or nrare 
environmental gradients Is called ordination (Goodall, 1954 cited In Kessell, 
1979). Since vegetation varies continuously along a moisture gradient, samples 
can Indeed be ordinated. 

Frequently, the development of a useful classification system requires the 
use of ordination methods. Discontinuities In the natural vegetation are 
sought for the purpose of determining the boundaries of the community types 
recognized. These are often best determined objectively by employing the 
methods of gradient analysis and ordination. The development of a Montana 
habitat-type system (Phlster et a1 ., 1977, cited In Kessell. 1979) Is a good 
example of the successful use of ordination In developing a classification 
system. 

Gradient modeling has been the first extensive application of gradient 
analysis to the needs of resource managmnent Information systems (Kessell. 
1979). Gradient modeling Involves tj-.? linkage of a multidimensional gradient 
analysis with a remote site-specific '/iventory and appropriate computer soft- 
ware. Once the gradient model Is CMiplete, It can provide quantitative com- 
munity Inferences (I.e.. biomass, cover) If the location of each site within 
the gradient matrix Is known (geographic coordinates, elevation, aspect, etc.}. 
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The Initial step Is to obtain Information about the vepetatlon. Data on 
the vegetation can be obtained by field samples ("ground truth" studies) 
and remote methods (aerial photography and satellite Imagery). Most systems 
use both. Detailed ground truth data are used to derive community- types, 
whereas aerial photographs and Imagery are generally used to Infer the 
vegetation of unsampled areas. 

To date, a considerable number of vegetation Investigations have been carried 
out using Landsat MSS Imagery but very little has been attempted with the 
meteorological satellite systems, particularly the NOAA/AVHRR. Gray and 
McCrary (1980, cited In Gray and McCrary, 1981) obtained a high correlation 
for detection of vegetation greenness between the NOAA-6 AVHRR Large Area 
Coverage (LAC) data sets and Landsat F6S data within Identical target areas. 

This finding led Gray and McCrary (1981) to suggest the NOAA satellite systems 
should be used for monitoring global vegetation. One major advantage of NOAA 
over Landsat Is the tremendous Increase In frequency of data collection. Gray 
and McCrary (1981) anticipated that variations In the AVHRR responses will 
provide Information about reactions of vegetation to moisture availability 
and thermal effects. They have demonstrated this for croplands In southern 
Texas before and after Hurricane Allen In April 1980. Since v.qetatlon, 
particularly In regions arid and semi-arid Is very responsive to moisture 
patterns. It Is well worthwhile to Investigate temporal changes In natural 
vegetation and how closely these relate to shifts In the AVHRP vegetation 
Index. 

Quite likely, we shall ultimately dscover that the success of stratifying 
different vegetation types from AVHRR vegetation Indites will depend not on 
spatial distinctions but on temporal distinctions, I.e., the rate and mag- 
nitude of the spectral shift during a single season. Eventually, vegetation 
Indices can be ordinated (Indirect gradient analysis) and correlated to 
ground truth vegetation and climatological graJlents (direct gradient analysis). 
The ultimate gradient model. Incorporating both field and satellite data, may 
permit vegetation classification and monitoring of changes with minimal ground 
truthlng. 
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Data Acquisition and Processing 

The vegetation gradient model Initially necessitated establishing a sample 
series along an enviroranental gradient. At NASA/JSC, we were geographically 
sitting at the eastern edge of perhaps one of the best natural east-west 
gradients In North America: Trans-Texas, along approximately 30 (°N) 
latitude. As one moves from Beaumont to El Paso, Texas, one passes through 
four major natural vegetation regions: (a) mixed pine-hardwood forest, 

(b) savannah, (c) shrubland, and (d) desert/desert scrub (Map 7). Pre- 
cipitation exhibits a continuun that has an annual mean >50" (east Texas) 
to<8" (west Texas). The elevation Is 0.0' at the Gulf, 5000' just east 
of El Paso. It would be difficult to Identify a better east-west continuum 
anywhere that changes gradually, yet dramatically and without any obvious 
dis junctures over a distance of approximately 750 miles. While the soil 
and geology definitely change across Texas, we do not Intend to include a 
discussion of those variables at this time. 

Our main objective has been to design a model that may ultimately allow 
vegetation classification on a global scale utilizing satellite imagery. 

We initially expected to use Landsat data. By a stroke of good fortune, 
we discovered that NOAA/AVHRR "Metsat" data was not only being archived 
locally by NOAA personnel (T. Gray, D. McCrary) in a readily useable form 
but it fit our specifications perfectly. The appropriate software had been 
written by Lockheed, Inc., to be able to retrieve raw pixel data for AVHRR - 
channels 1 and 2 along specific scanlines or bands of scanlines across the 
entire state of Texas. In addition, the software provided geographic 
coordinates for each pixel. In order to obtain a specific scanline, it 
was really only necessary to provide the specific coordinates at the be- 
ginning and end of our trans-Texas transect. Since the NOAA - n series 
of satellites orbit is near-polar, sun-synchronous, and twice daily, it 
crosses a given longitude at an angle and at varying places. Consequently, 
since scanlines are perpendicular to the orbit, it was impossible to select 
scanlines that remained "isolatitude" or were exactly superimposed from one 
date to the next. It Is also important to remember that NOAA scenes cover 
such an expanse that related angles to each pixel vary greatly. To permit 
comparisons, the pixel radiance values have been normalized to an overhead 
sun. 
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It has been our Intention to acquire scanllnes for four cloud- free days 
during 1980 • one from each season. At this point we have only been able 
to process 3 dates (April 19, July 10, October 9); winter has been excluded 
for lack of data. We requested and received bands of 5 adjacent scanllnes, 
extending essentially from El Paso to Beaumont (Figure 2). At predetermined 
locations which corresponded to our ground truth sites along the strip, we 
sampled a 25 pixel grid (5 x 5), obtaining an average grid value of pixel 
counts for each two channels. The selection of the 25 pixel sample grids 
was somewhat difficult because It was not possible to accurately ground- 
truth the transect line. The Intention has been to select 12 sites, approx- 
imately 3 sites In each of the four major vegetation regions bisected (Map 7), 
Using Texas vegetation-type maps (Texas Parks and 'Midlife Department - based 
on Landsat data), original Landsat MSG scenes, aerial photos. Aeronautical 
Navigation maps (1:1,000,000), and selected vegetation references (Gould, 

1975; see Smeins, 1978), an effort was made to choose "homogeneous" natural 
vegetation sites, devoid of water, urbanization, and cultivation. The site 
locations were shifted slightly between sampling dates because the scanllnes 
could not be superimposed. 

It was difficult deciding just how to Initially treat the satellite data. 

Gray and McCrary (1981) have devised their own vegetation index, that they 
now rather appropriately call the Gray-McCrary Index (GMI). The GMI is 
simply the difference between the solar-zenith corrected albedo value for 
the two channels. At least initially, we are using the Landsat-derived 
normalized difference (ND) equation of Rouse, et al. (1973) and Deering, 
et al. (1975) where: 

Channel 2 - Channel 1 
ND = 

Channel 2 + Channel 1 
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Figure 2.- Example of a band of 5 NOAA/AVHRR scanlines and the 
5x5 sampling grids along the Texas transect. 





Normalized difference values were obtained for each of the twelve (12) 
sites on the three (3) dates. The next step was to decide how to best 
statistically treat the ND values. At this stage we have simply been 
able to do some preliminary analysis of variance. Normalized difference 
values were plotted as a function of longitude, sponge and certain vegetation 
characteristics (i.e., biomass, net productivity, leaf area). 

Vegetation Regions on the Texas Transect 

The Texas transect selected for our model essentially runs from Beaumont 
(94°W) to El Paso (106°) and the 12 sample sites have been numbered east 
to west (see Map 7). Since it was not feasible to visit the transect for 
optimal site selection, it was necessary to utilize vegetation maps. At 
this level of our investigation, the ground-truth precision was not terribly 
critical since our initial concern has been tc get a general feeling for the 
potential of NOAA imagery for global vegetation stratification. 

Gould's (1975) vegetation map and discussion of the vegetation regions of 
the state is large-scale but is the best complete map available. The 
Texas Parks and Wildlife Department is in the process of completing a state- 
wide series of land-use classification maps based on Landsat data. Utilizing 
primarily those two sources, it appears as though our transect bisects four 
major vegetation regions. Table III enumerates those four regions from east 
to west, their approximate longitudinal boundaries on the transect, and 
document vegetation. For a more complete vegetation description, see Gould 
(1975), Texas Parks and Wildlife vegetation-type maps, and Smeins (1978). 
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TABLE III. THE VEGETATION REGIONS. THEIR APPROXIMATE LONGITUDINAL BOUNDARIES, AND 
DOMINANT VEGETATION ALONG THE TEXAS TRANSECT. 



f ** E 

•o 

•a 


32 

</) 5 

c 

• c 


o 

^ 10 Cft 

* 10 

k- <0 


*— » 

I— o.^ 

E •- 

0) f- 

.C *3 

h" 

</) u 

1— o. 

O. 0) 

«0 C 

< 

* <o 

V p 

*r“ 0> 

3 <0 



s. 

c c 

.O f— 

UJ 

(L* n 

•> o 

3 (0 

k. 0> 

C9 

C C7) 

W 

•o k- 

10 Ol 

Ui 


lO 


•O C 


a. Oi c 

o • 

A 

<0 


0) 10 

>» 

•C 

•> U 

»— 

>>5 0 

• &. 

iA 0) 

0) w- 

z 

r~ W 0) 

(/) C- 

<0 ■♦-> 


<c 

1 — a. 

o <u 


O • 


o • 

Kn .Q 

• 3 

Vi <0 

HH 

(U * 

V) 

Vi O' 

O U 


^ C E 

<0 o 

10 

0 ) o 

o 

O •«- 1— 

U (O 

10 <u 

k. 3 

o 

r- Ol OI 

o»jc 

O E 

O >> 


o 

o 




LO 

o 



o 

•O’ 

CO 

o 



• 

• 

CM 


UI (/) 

LO 

00 

o 

O 

Q h- 

cr> 

cr» 


VO 


1 

i 

1 

o 

h- s: 

o 

o 

o 


•—1 *— t 

in 

in 

o 

1 

O _l 



CO 

o 


• 

• 

• 

CM 

O 

ro 

in 

00 

o 


<r> 

a> 

<y» 




•» m 


to 



x: OJ 


c 

>ui 


(O 



m o£. 


c u 

3 

Vi 

r>*. r> 


C *r- 

<o 

(O 

h- 


<0 (0 


m 

r— < 


> U 

-4-> 


^ _J 


lO Q. 

<o 

to eO 

O 

t/> 

(/) 


o 

o z 

•a 

■o 

CL 

U to 

-J UI 

o 

c 


Oi c 

=3 z: 

o 

lO <o 

V) 

U. *p- 

o o 

5 

O f— 

•o 

<0 

CD Z 

>> 


u 

to -M 



-»-> o 

(0 

C C 


c 

10 «0 

5 

«0 3 


'T— 

O f— 

-o 



Q. 

0_ CQ 

UI 

H- it 






o 





H- « 





CD 




1 

UI 

■o 



+J 

OC 

O 



k. 


O 



0) 

z 

5 



Vi 

o 

T3 


•o 

0) 


U 4-> 

JO. 

c 

a 

»— 

<0 (/) 

ID 

(0 

— ..Q 

< 

z <u 

c 

r— 

4-> 3 


1 s- 

c 

.n 

k. k. 

UI 

<u o 

lO 

3 

<u a 

CD 

C U- 

> 

U 

to to 

UI 


<o 

JC 

0) 

> 

Cl 

CO 

CO 

o 


26 






i 


Results and Discussion 


1. Norwelized difference (ND) as a function of longitude . As pre- 
viously mentioned It was Impossible to Insure superimposition of pixel 
sampling locations between dates because the satellite orbit fluctuates. 

For example, while still In pine- hardwood forest, site 1 on April 19 Is 
not geographically Identical to site 1 on October 9. The tabular and 
graphic summaries when the three sample dates were Individually plotted 
as a function of longitude are presented below (Table IV, Figures 3, 4). 

There Is generally a high correlation between ND and longitude 
(mean r = .756); NO decreases from east to west. 

The sample size of 12 Is not sufficiently large enough to merit serious 
discussion as to significant differences among the three sampling periods. 
Neither can we eliminate the real possibility of cloud cover affecting 
reflectance values. We definitely know that there was some cloud cover 
In West Texas on October 9. What the regressions do suggest Is what one 
would expect knowing the phenological nature of the vegetation regions on 
the transect. April and July values are high In the mixed forest because 
the deciduous trees have leafed out. By October, they have dropped their 
leaves, reducing their "greenness”. Progressing westward across the state, 
a greater percentage of the perennials are non-declduous but the vegetation 
becomes less dense and more dependent on Infrequent precipitation. In 
addition, the amount of exposed ground In the desert scrub poses problems 
of separating soil spectra from vegetation spectroreflectance (Miller, 

Lee D. , pers. coiwn,). 

The variation between seasons appears to diminish going from east to west. 
The mean regression for all three dates (Figure 4) Indicates not only that 
the greater between-date variation Is In the pine-hardwood forest but also 
that April, and to a lesser extent July are well above the line for the 
three dates. We would have expected much higher values 'Juring the annual 
blooms In April In the Trans-Pecos but more extensive sampling may clarify 
this. 
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TABLE IV. INDIVIDUAL NORMALIZED DIFFERENCE (ND) VALUES. MEANS, REGRESSION 
EQUATIONS, AND FOR 3 DATES IN 1980 ALONG TEXAS TRANSECT. 


Site April 19 


1. 

.258 

2. 

.289 

3. 

.286 

4. 

.259 

5. 

.224 

6. 

.126 

7. 

.138 

8. 

.062 

9. 

.038 

10. 

.044 

11. 

.046 

12. 

.035 


/ 


July 10 October 9 


238 

.161 

229 

.171 

231 

.215 

136 

.147 

205 

.151 

137 

.090 

125 

.033 

141 

.101 

159 

.112 

029 

.051 

005 

.044 

001 

.063 


Combined 


A 

Y 


.159 

2.526 - .024X 
.752 


.128 

2.268 - .022X 
.934 


.112 

1.23 - .OllX 
.636 


.130 

2.00 - .019X 
.756 
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Figure 4.- 1980 Vegetation indices and mean regression line - Texas transect. 



2. Sponge Index as a function of Ibngltudt . For a discussion of a 
iiw moisture variablet the sponge, we refer you to earlier sections II, III 
of this report. Long-term annual “normals" (1941-70) were calculated for 
26 stations along or near the Texas transect. These values were plotted 

as a function of longitude (Figure 5). The results Indicate an extremely 
high correlation (r ■ .911) between sponge and longitude; sponge increases 
(as does precipitation) from west to east. 

3. Normalized difference (ND) as function of sponge . Since both the 
vegetation index, ND, and sponge showed similar positive correlation with 
geographic position on the transect, it seemed appropriate to interpolate 
sponge values from the 26 stations (see Table A**, page 54; Maps 8,9) 
along the transect for each of 36 ND values (3 values for each of the 12 
sites). When ND was plotted against sponge, there was a good correlation 
(r “ .777) (Figure 6). Regressing the two regressions, (ND (Figure 4) 

vs. Sponge (Figure 5), the result is a ' ery pleasing (r^- 1.001) but transect- 

limited, prediction model that permits estimating longitude, sponge index, 

and ND value, requiring input of only one of the 3 variables (Figure 7). 

This model serves to illustrate the very high correlation between the 

vegetation index and the sponge. To further establish this correlation, 

if one plots the highest ND value of the three dates at each 12 stations 

2 

against the long-term sponge value for that particular month, r ■ .946. 

4. Normalized difference (ND) as a function of vegetation . As earlier 
mentioned, it was not feasible to actually ground truth the transect to 
verify and quantify the vegetation. No doubt this should be done at some 
later stage. In lieu of a better alternative, biomass, net productivity, 
and leaf area estimates from Whittaker and Likens (1973) were used (see 
Table V). Since the predominant vegetation limiting factor along the Texas 
gradient is moisture, these vegetation parameters predictably decline from 
east to west. It has been previously demonstrated that all of these para- 
meters have been correlated with spectral data (see introduction to Tucker, 
et al. 1981 for literature review). While general, when ND means for each 
site are plotted against biomass, net productivity, and leuf area mean 
values, the results are prophetic as to which vegetation characteristic has 
the highest correlation with the satellite data (Table VI). 
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TABLE V. MEAN VALUES FOR WORLD-WIDE ESTIMATES OF NET PRODUCTIVITY BIOMASS 
AND LEAF AREA INDEX (WHITTAKER AND LIKENS. 1973) 



NET PRODUCTIVITY 
om/M^/Yr 

BIOMASS 

Kg/M^ 

LAI 

Temperate Evergreen/Deciduous 
Forest 

1250 

32.5 

8.5 

Savannah 

900 

4.0 

4.0 

Shrubland 

700 

6.0 

4.0 

Desert/Desert 

90 

0.7 

1.0 


32 


TABLE VI.- COEFFICIENTS OF DETERMINATION (r^) WHEN MEAN ND VALUES FOR THE 12 
TEXAS TRANSECT SITES ARE PLOTTED AGAINST NET PRODUCTIVITY, LEAF 

AREA. AND BIOMASS 


in 

CSi Cft 

1 - 00 


00 


o 

in 

iO 
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Figure 7.- Vegetation-Sponge prediction model - Texas transect. 
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HAP 8. STATION LOCATIONS FOR SPONGE CALCULATIONS ALONG TEXAS TRANSECT 
(PERIOD 1941-70 AND 1980, EXCEPT AS NOTED IN TABLES AI (a) - (f)) 
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MAP 9. MEAN 8- INCH SPONGE VALUES ALONG TEXAS TRANSECT 7 

(1980) 
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The Inference is that net productivity has the highest correlation to the 
vegetation index, biomass the lowest. Several studies have recently shown 
that currently used remote sensing techniques are not sensitive to non-green- 
leaf components of the phytomass. However, there appears to be a high corre- 
lation of spectral data with green-leaf area (biomass) and net production of 
certain vegetation types (see Introduction, Tucker, et al. 1981). Our selection 
of NO as a vegetation index was largely founded on Oeering and Haas' (1977) high 
correlation between Landsat-derived ND and rangeland biomass. 

In Table V, it is noted that of the three parameters of vegetation, biomass 
is the only one that doesn't consistently decline from east to west along the 
transect. Shrubland, consisting largely of woody perennials, does not produce 
the annual net production that a savannah, containing more herbaceous annuals 
would, but its accumulative biomass would be greater. The ND does not re- 
spond to the increase in biomass from savannah to shrubland because much of 
that biomass is tied up in non-green components in shrubland which is not as 
true in the savannah. Because net productivity and leaf area more closely re- 
flect the actual spectral component of those vegetation regions on our transect, 
we would anticipate their higher correlation with ND. 

5. Vegetation-Sponge Index (VSI) . Having previously established a high 
correlation between ND and sponge, we would like to propose a new index that 
represents the multiplicative of the two variables: the Vegetation - Sponge 

Index (VSI). 


ND X SPONGE = VSI 

The mean VSI values for the 12 sample sites are presented below in tabular 
and graphic form (Table VII, Figure 8). 
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TABLE VII. SPONGE, ND, AND VSI MEAN VALUES FOR DATES IN 1980 (APRIL 19. JULY 10. 
OCTOBER 9) FOR 12 SITES ALONG TEXAS TRANSECT. 
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ND - NORMALIZED DIFFERENCE 








When VSI Is plotted as a function of longitude, the separation of the four 
major vegetation regions becomes much more apparent (Figure 8). It Is 
ever possible to suggest boundaries for the four classes: 

’ >1.0 pine-hardwood forest 

0. 4-1.0 savannah/cropland 

.05-0.4 shrubland 

< .05 desert/desert scrub 

Using these suggested limits. It becomes readily apparent that the least 
within-class variation occurs In desert and forest; the greatest variation 
is in the shrubland and savannah/cropland classes. This seems consistent 
with what we would expect. The desert and mixed forest would be more homo- 
geneous In the sense that desert has extensive bare soil during much of the 
year while the forest would have relatively little. The shrubland and 
savannah would be considerably more heterogeneous. 
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V. THE TEXAS MODEL; CONCLUSIONS AND PROJECTIONS 


It has batn dcmonstritad that the sponge variable Is a superior tool for 
the analysis of cllmate/vegetatlon relationships. Furthermore, NOAA/AVHRR 
satellite data proved useful for vegetation stratification. Finally, a 
preliminary multivariate model of vegetation distribution, the VSI, was 
developed based on an experimental east-west Texas gradient. The next 
stage of this research effort will Involve a more extensive analysis of 
the application of the model to the Texas transect (e.g.. Increasing the 
sample size), to be followed later by refinement of the model which will 
then be tested against other natural vegetation regions across North 
America. It Is anticipated that ultimately such a model may be utilized 
for global vegetation surveys, and as a means of remotely monitoring 
vegetation region dynamics (e.g., desertification). 
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Except as noted. All values in Inches. Data from Climatological Summaries for Texas. National Heather Service. NOAA. 
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Except «s noted. A\1 values In Inches. Data from CUmatolo firal Sunmarles for Texas, National Heather Service 
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Data froai Cllawtoloqlcai Suaaaarles for Texas, «tat1onal Weather Service. 
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Eactpt «t noted. All valuts In Inches. Data fn»a C11«iatoloq1ca1 Stanaarlas for Texas. National Weather Service. 




































































































































































































































TABLE Al(f). PERIODS OF RECORD OTHER THAN 1941-1970: TABLES 2(a) - ’(e) 


’l931 - 1969 
^1931 - 1970 
^932 - 1969 
^933 - 1967 
^1933 - 1968 
^1934 - 1969 
^937 - 1966 
®1937 - 1969 
®1938 - 1967 
’®1939 - 1968 
^^939 - 1969 
’^1940 - 1966 
^®1940 - 1969 
^^1941 - 1967 
^®1941 - 1968 
^®1942 - 1971 
^^943 - 1970 


’®1947 - 1966 
’®1947 - 1970 
^®1948 - 1967 
^^1948 - 1969 
^^1948 - 1970 
^®1949 - 1966 
^^1949 - 1967 
^®1949 - 1968 
^®1949 - 1969 
^^949 - 1970 
^®1950 - 1969 
^®1950 - 1971 
^°1951 - 1966 
®’l951 - 1967 
®^1953 - 1970 
®^1956 - 1970 
®'*1959 - 1971 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


TABLE All. PRECIPITATION AND 8-INCH SPONGE VALUES AT SELECTED STATIONS 
ALONG AN EAST-WEST TEXAS TRANSECT; BEAUMONT (30.08ON; BA.IO^W) 


Station 


Sta. 

No. 

Period 

Lati- 
tude (°N) 

Longi- 
tude (''Ul 

Januai 

PrecTpTH 

^nge 

July 

PreclD. Sponge 

Year 

Precip. 

Sponge 

Beaumont 


0613 

1941-70 

30.08 

94.10 

4.57 

7.00 

5.71 

4.30 

54.77 

5.20 

Beaumont 


0613 

1980 

30.08 

94.10 

5.84 

5.04 

0.95 

1.34 

57.59 

4.72 

Beaumont 


0613 

1979 

30.08 

94.10 

7.91 

7.43 

15.50 

3.93 

79.56 

5.55 

Lufkin 


4524 

1941-70 

31.23 

94.75 

4.29 

6.32 

2.83 

2.63 

44.93 

4.22 

Liberty 


5196 

1941-70 

30.05 

94.82 

4.21 

6.21 

4.65 

3.49 

49.61 

4.41 

Liberty 


5196 

1980 

30.05 

94.82 

0.40 

3.99 

0.36 

2.65 

56.77 

4.84 

Liberty 


5196 

1979 

30.05 

94.82 

8.91 

7.43 

7.61 

4.08 

70.11 

5.31 

Huntsville 


4382 

1941-70 

30.72 

95.57 

3.80 

6.25 

3.28 

3.10 

45.95 

4.39 

Huntsville 


4382 

1980 

30.72 

95.57 

0.94 

2.61 

0.05 

0.78 

28.30 

2.67 

Huntsville 


4382 

1979 

30.72 

95.57 

7.65 

7.64 

5.90 

3.20 

64.07 

5.31 

Brenham 


1048 

1941-70 

30.15 

96.40 

2.78 

4.34 

1.90 

2.01 

38.96 

3.53 

Temple 


8910 

1941-70 

31. in 

97.35 

2.35 

4.08 

1.96 

1.97 

33.87 

3.03 

San Marcos 


7983 

1941-70 

29.88 

97.95 

2.06 

3.08 

1.89 

2.04 

33.86 

2.74 

New Braunfels 

0832 

1941-70 

30.10 

98.42 

2.12 

3.23 

1.98 

1.79 

34.39 

2.86 

Blanco 


6276 

1941-70 

29.70 

98.12 

1.88 

2.92 

1.83 

1.78 

32.61 

2.56 

Llano 


5272 

1941-70 

.30.75 

98.68 

1.37 

2.16 

1.20 

1.22 

26.16 

2.12 

Brownwood 


1138 

1941-70 

31.72 

98.98 

1.72 

2.63 

1.85 

1.92 

27.20 

2.32 

Korrvllle 


4782 

1941-70 

30.05 

99.15 

1.86 

2.36 

2.10 

1.96 

30.75 

2.59 

Brady 


1017 

1941-70 

31.12 

99.35 

1.52 

2.16 

1.34 

1.41 

23.27 

1.96 

Brady 


1017 

1980 

31.12 

99.35 

1.18 

2.42 

0.00 

0.44 

24.21 

2.11 

Brady 


1017 

1979 

31.12 

99.35 

1.14 

2.98 

1.96 

1.33 

23.47 

2.09 

Sonora 


8449 

1941-70 

30.57 

100.65 

0.82 

1.C8 

1.61 

1.33 

19.28 

1.46 

Colorado City 

4974 

1941-70 

32.38 

100.87 

0.82 

1.46 

2.14 

1.39 

19.80 

1.44 

Otona 


6734 

1941-70 

30.72 

101.20 

0.81 

1.21 

1.30 

1.19 

17.59 

1.36 

Ozona 


6734 

1980 

30.72 

101.20 

0.08 

2.52 

0.00 

0.87 

17.10 

1.61 

Bio Spring 


0786 

1941-70 

32.25 

101.45 

0.63 

1.30 

1.97 

1.28 

15.72 

1.28 

McCamey 


5707 

1941-70 

31.13 

102.20 

0.64 

0.98 

1.64 

1.42 

12.75 

0.94 

He Carney 


5707 

1979 

31.13 

102.20 

0.08 

0.67 

1.01 

0.58 


0.57 

Fort Stockton 

3278 

1941-70 

30.87 

102.92 

0.83 

0.98 

1.44 

0.95 

12.23 

0.86 

Mink 


9829 

1941-70 

31.78 

103.20 

0.59 

0.79 

1.64 

0.86 

11.11 

0.79 

Pecos 


6892 

1941-70 

31.42 

103.50 

0.34 

0.57 

1.35 

0.66 

9.05 

0.58 

Pecos 


6892 

1980 

31.42 

103.50 

0.04 

1.46 

0.09 

1.44 

15.78 

1.12 

Pecos 


6892 

1979 

31.42 

103.50 

0.75 

1.10 

0.94 

0.40 

8.19 

0.57 

Balmorhea 


0498 

1941-70 

30.98 

103.75 

0.64 

0.86 

1.56 

0.86 

12.12 

0.81 

Balnorhea 


0498 

1980 

30.98 

103.75 

0.26 

1.27 

0.00 

0.14 

13.15 

1.45 

Balmorhea 


0498 

1979 

30.98 

103.75 

1.43 

1.51 

1.89 

0.68 

12.78 

0.87 

Van Horn 


9311 

1941-70 

31.05 

104.83 

0.54 

0.89 

1.75 

0.80 

10.23 

0.78 

Salt Flat 


7920 

1941-70 

31.78 

104.90 

0.30 

0.65 

1.48 

0.74 

8.51 

0.66 

Salt Flat 


7920 

1980 

31.78 

104.90 

0.11 

0.56 

0.16 

0.10 

7.08 

0.56 

Salt Flat 


7920 

1979 

31 ./u 

104. DO 

0.88 

3.68 

2.65 

0.79 

9.47 

0.65 

Ysleta 


9966 

1941-70 

31.70 

106.32 

0.40 

0.56 

1.64 

0.86 

11.11 

0.49 

Ysleta 


9966 

1930 

31.70 

106.32 

0.96 

1.38 

0.00 

0.03 

9.30 

0.38 

Ysleta 


9966 

1979 

31.70 

106.32 

0.80 

0.51 

1.37 

0.49 

8.17 

0.52 

La Tuna 


4931 

1941-70 

31.97 

106.60 

0.41 

0.71 

1.61 

0.66 

8.C‘: 

0.59 

'a 11 oreclc 
Ashevi I 

iltatlon 
1e, N.C. 

and sponge values In 

Inches. 

Data from National 

Climatic 

Center, 

NOAA, 
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Texts Sponge Values: 19:io’ 


West 



I <JSt 



'oased on Z8 stations along l-W “Texas Transect" from lieaunont to La Tuna. 
1941-1970; 7*j Stations throughout Texas. 

.'lote; ?5; Vertical Exaggeration on all plots. 


Figure Al.- Texas sponge values for 1980 compared to 
normal annual Texas sponge values. 
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Viilu<‘'<: Aiifll. I'HKI 



llormiil Icx.i'. '.IMIII'H’ V.iliic ; li'ril 



on '/X '.(.il ioir. .iIc«h| I -W “Irx.i'. 

^ I'I'1 1 - 1')/0; /'■ '.liii iiin'. Ihniui|hciiil, li'x.i 
Noll-; '/'<> Vcrlii.il i'x.iiii|cr.ilii)n on .ill 


Ir.iiiM’i I.” IriBii Ili'.iimxiMl lol.i lin.i. 


plol',. 


Figure 42.- Texas sponge values for April 1980 compared to 
normal Texas sponge values for April. 
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ORIGINAL PAGE 
OF POOR QUAUTY 


Irx.is Vtilurs; .Inly, I'lJld 



tkinatl ti‘X<is S|HiniK' V<ilui-s: 





*U.i«ird on ^It Millions <iloni| I -II “Tcxiis Trtiiiscil" froNi IkMimonl to I o I'ln.. 
''l9d l-t'l/O; ?!i slot ions throiKihoul Toxos. 

Note'; X!r.' Vert i (111 I Xiii|(|criilion on iill |ilols. 


Figure A3.- Texas sponge values for July 1980 compared to 
normal Texas sponge values for July. 
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ORIQ'.NAL r Aui: ' 
OF POOR QUALITY 


Texts Sponge Vtiucs: October, 19)10* 



Nomdl Texts Sponge Vtlues: 


October^ 


West 


Fast 


Based on 28 stations along F-M "Texts Transect" froxi Beaunont to La Tuna. 
^)9t1-)970 ; 7S stations t)iroug)iout Texas. 


itote: ?5T. Vertical Exaggeration on all plots. 



Figure A4.- Texas sponge values for October 1980 compared 
to normal Texas sponge values for October. 
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